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ABSTRACT 

Recently, many dark energy models whose equation-of-state parameter can cross the phantom 
divide Wde = —1 have been proposed. In a previous paper [Class. Quant. Grav. 22, 3189 (2005); 
hep-th/0501160 , we suggest such a model named hessence, in which a non-canonical complex scalar 
field plays the role of dark energy. In this work, the cosmological evolution of the hessence dark 
energy is investigated. We consider two cases: one is the hessnece field with an exponential potential, 
and the other is with a (inverse) power law potential. We separately investigate the dynamical system 
with four different interaction forms between hessence and background perfect fluid. It is found that 
the big rip never appears in the hessence model, even in the most general case, beyond particular 
potentials and interaction forms. 
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I. INTRODUCTION 

Dark energy Q] has been one of the focuses in modern cosmology since the discovery of accelerated 
expansion of the universe 0, IE 0, IE IE HE| ■ The simplest candidate of dark energy is a tiny positive 
cosmological constant. As an alternative to the cosmological constant, some dynamical scalar field models 
have been proposed, such as quintessence 0, E3> phantom [E |E 0, k-essence etc. In the 
observational cosmology of dark energy, equation-of-state parameter (EoS) Wde = Pde/ Pde plays a central 
role, where pd e and pd e arc the pressure and energy density of dark energy respectively. The most 
important difference between cosmological constant and dynamical scalar fields is that the EoS of the 
former is always a constant, —1, while the EoS of the latter can be variable during the evolution of the 
universe. 

Recently, by fitting the SNe la data, marginal evidence for Wd e (z) < — 1 at z < 0.2 has been found [l2T |. 
In addition, many best-fits of the present value of Wde are less than —1 in various data fittings with 
different parameterizations (see |l3| for a recent review). The present data seem to slightly favor an 
evolving dark energy with Wde being below — 1 around present epoch from Wde > — 1 in the near past . 
Obviously, the EoS Wd, F cannot cross the so-call ed p hantom divide Wde . = — 1 for quintessence or phantom 
alone. Some efforts [3 13 13 13 13 EE EE El EE El EE EE Ell El have been made to build dark 
energy model whose EoS can cross the phantom divide. Here we mention some of them, such as the 
geometric approach (brane or string) [l8|. holographic dark energy 0, l2C| . and scalar-tensor theory 
(esp. non- minimal coupled scalar field to gravity) (221 etc. 

Although some variants of k-essence look possible to give a promising solution to cross the phantom 
divide, a no-go theorem, shown in |2l|, shatters this kind of hopes: it is impossible to cross the phantom 
divide Wde = — 1) provided that the following conditions are satisfied: (i) classical level, (ii) general 
relativity is valid, (iii) single real scalar field, (iv) arbitrary Lagrangian density p(ip,X), where X = 
ig^dpcpdvip is the kinetic energy term, (v) p(tp, X) is continuous function and is differentiable enough, 
and (vi) the scalar field is minimal coupled to gravity. Thus, to implement the transition from Wd e > — 1 
to Wde < — 1 or vice versa, it is necessary to give up at least one of conditions mentioned above. 

Obviously, the simplest way to get around this no-go theorem is to consider a model with two real 
scalar fields, i.e. to break the third condition. Feng, Wang and Zhang in [l4| proposed a so-called quintom 
model which is a hybrid of quintessence and phantom (so the name quintom) . Naively, one may consider 
a Lagrangian density [3, El EE El 

1 2 l 2 

£ quintom = -^{d^i) — - {8^ fa) ~ V {fa , (fa), (1) 

where fa and fa are two real scalar fields and play the roles of quintessence and phantom respectively. 
Considering a spatially flat Friedmann-Robertson- Walker (FRW) universe and assuming the scalar fields 
fa and fa are homogeneous, one obtains the effective pressure and energy density for the quintom 

Pquintom = ~ \^>\ ~ V(fa, fa ), Pquintom = ]z<t>\ " \^>\ + V(fa,fa), (2) 

respectively. The corresponding effective EoS is given by 

_ # -4>l- 2V(fa, fa) 

Wquintom -^-4 + 2V(fa,fay {d) 

It is easy to see that w qu i n tom > — 1 when 4>\ > 4>\ while w qu i n tom < — 1 when <j>\ < 4>2- The cosmological 
evolutions of the quintom model without direct coupling between fa and fa, i.e. V(fa,fa) = V ( f >1 + 
Vfc = V^jO exp(— XfaKfa) + V$ 2 o exp(— X,p 2 Kfa), and with a special interaction between fa and fa, i.e. 

V(fa,fa) = Vfa + Vfa + Vi n t and Vi nt ~ (Vfa V^) 1 ^ 2 , were studied by Guo et al and Zhang et al [3], 
respectively. They showed that the transition from w qu i n tom > — 1 to w qu i ntom < — 1 or vice versa is 
possible in this type of quintom model. 

In |28[ , by a new view of quintom dark energy, we proposed a novel non-canonical complex scalar field, 
which was named "hessence" , to play the role of quintom. In the hessence model, the phantom-like role 
is played by the so-called internal motion 9, where 8 is the internal degree of freedom of hessence. The 
transition from wt > —1 to Wh < — 1 or vice versa is also possible in the hessence model |28j. We will 
briefly present the main points of hessence model in Sec. [n] 
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The main aim of this work is to investigate the cosmological evolution of hessence dark energy. We 
consider the hessence and background perfect fluid as a dynamical system [2^. By the phase-space 
analysis, we find that some stable attractors can exist, which are either scaling solutions or hessence- 
dominated solutions with EoS (wh and w e ff) larger than or equal to —1. No phantom- like late time 
attractors with EoS less than —1 exist. This result is very different from the quintom model considered 
in |22l |23|. where the phantom-dominated solution is the unique late time attractor. If the universe is 
attracted into the unique phantom-dominated attractor with EoS less than —1, it will un derg o a super- 
accelerated expansion (i.e. H > 0) forever and the big rip is inevitable (see also HHaall)- On 
the contrary, in the hessence model, EoS less than —1 is transient. Eventually, it will come back to 
quintessence-like attractors whose EoS is larger than —1 or asymptotically to de Sitter attractor whose 
EoS is a constant — 1. Therefore, the big rip will not appear in the hessence model. 

The plan of this paper is as follows. In Sec.[n] we will briefly present the main points of the hessence 
model. In Sec. lIIII we give out the equations of the dynamical system of hessence with/without interaction 
to background perfect fluid for the most general case. That is, we leave the potential of hessence and the 
interaction form undetermined. We will investigate the dynamical system for the models with exponential 
and (inverse) power law potentials in Sec. lIVI and Sec.0 respectively. In each case with different potential, 
we consider four different interaction forms between hessence and background perfect fluid. The first one 
corresponds to the case without interaction, i.e. the interaction term is zero. The other three forms are 
taken to be the most familiar interaction ones considered in the literature. In all these cases, we find 
that no phantom- like late time attractors with EoS less than —1 exist. In Sec. IVII we will show the big 
rip will not appear in a general case beyond the above considered cases. Finally, brief conclusion and 
discussions are given in Sec. IVIII 

We use the units % — c = 1, n 2 = 8irG and adopt the metric convention as (+, — , — , — ) throughout 
this paper. 



II. HESSENCE DARK ENERGY 

Following [2^, we consider a non-canonical complex scalar field as the dark energy, namely hessence, 

$ = 0i+i02, (4) 

with a Lagrangian density 

c h = \[ {d^f + (d,<£*) 2 ] u(<p 2 + <r 2 ) = \ [ {d^f tfid.ef } - v(<t>), (5) 

where we have introduced two new variables (0, 9) to describe the hessence, i.e. 

0i = 0cosh0, 2 = 0sinh6>, (6) 

which are defined by 

cf> 2 = <pj - <f>l, coth0=^-. (7) 

02 

Considering a spatially flat FRW universe with scale factor a(t) and assuming and 6 are homogeneous, 
from Eq. (jSJ) we obtain the equations of motion for and (9, 

4> + 3H<j> + 4>9 2 + V 4> = 0, (8) 
<f) 2 6 + (2(f»j) + 3H<t> 2 )6 = 0, (9) 

where H = a/a is the Hubble parameter, a dot and the subscript denote the derivatives with respect 
to cosmic time t and 0, respectively. The pressure and energy density of the hessence are 



PH 2 



1 (0 2 - 4?e 2 ) - Ph = \ (> - 2 ^ 2 ) + v(<t>), (io) 

respectively. Eq. © implies 

Q = a 3 <j} 2 9 = const. (11) 
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which is associated with the total conserved charge within the physical volume due to the internal 
symmetry j2^|. It turns out 



g_ Q 

a 3 (j) 2 

Substituting this into Eqs. (JBJ) and l)10[l. we can recast them as 



(12) 



+ 3H^+-^ + V rk = 0, (13) 



It is worth noting that Eq. (|13f) is equivalent to the energy conservation equation of hessence ph + 
3H (ph + Ph) = 0. The Friedmann equation and Raychaudhuri equation are given by, respectively, 

H 2 = y (Ph + Pm) , (15) 



H = ~—(ph + Pm+Ph+Pm) , (16) 

where p m and p m are the pressure and energy density of background matter, respectively. The EoS of 
hessence Wh = Ph/ Ph- It is easy to see that Wh > —1 when (j) 2 > Q 2 /(a 6 (j) 2 ) while Wh < —1 when 
tj) 2 < Q 2 /(a 6 (f) 2 ). The transition occurs when <p 2 = Q 2 /(a Q cj> 2 ). 

There are other interesting points in the hessence model, such as the avoidance of Q-ball formation, 
the novel possibility of anti-dark energy, the possible relation between hessence and Chaplygin gas, etc. 
We refer to the original paper [2^ for more details. 



III. DYNAMICAL SYSTEM OF HESSENCE WITH/WITHOUT INTERACTION TO 

BACKGROUND PERFECT FLUID 

Now, we generalize the original hessence model t ,28] to a more extensive case. We consider a universe 
containing both hessence dark energy and background matter. The background matter is described by a 
perfect fluid with barotropic equation of state 

Pm = W m p m = (7 - l)p m , (17) 

where the so-called barotropic index 7 is a constant and satisfies < 7 < 2. In particular, 7 = 1 and 4/3 
correspond to dust matter and radiation, respectively. 

We assume the hessence and background matter interact through an interaction term C, according to 

Ph + 3H(p h + Ph ) = -C, (18) 

p m +3H (Pm+Pm) = C, (19) 

which preserves the total energy conservation equation p to t + 3iJ (ptot + Ptot) = 0. Clearly, C = 
corresponds to no interaction between hessence and background matter. It is worth noting that Eq. I|13|) 
should be changed when C ^ 0, a new term due to C will appear in the right hand side. Since 9 is 
the internal degree of freedom [28], Eqs. ©, i|ll|l and 1)12(1 still hold. The hessence interacts to external 
matter only thr ough <j>. Thus, the equation of motion of <j> should be changed by the interaction. 
Following |3lL l32l l33| , we introduce following dimensionless variables 

_ K(f) _ kW _ K^fp~m~ _ V6 _ K Q 
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By the help of Eqs. I|14|) (|16[) . the evolution equations l|18fl and l|19|) can then be rewritten as a dynamical 
system: 



x' = 3x(x 2 -v 2 + j-z 2 -l)-uv 2 - -^A- - d, (21) 



7 2 A 2 kV ,4 
—z — 1 1 — uv =— 

2 > 

y = 3y [x —v + —z H ?= ^tt, (22 

V 2 / 2V3HVVH 7 v y 

^' - 3z (x 2 -v 2 + lz 2 ~ Tj + C 2 , (23) 

2 



a;tt , (24) 

7, 
2' 



3« r 



x 2 - v 2 + -z 2 -1) - xuv, (25) 



where 



C = K ° C ee K ° (26) 

a prime denotes derivative with respect to the so-called e-folding time jV= In a, and we have used 

11 o( x 2_ v 2 + 7 2 y (27) 



i/ 2 V 2' 
The Friedmann constraint equation H15J) becomes 

x 2 + y 2 + z 2 - v 2 = 1. (28) 
The fractional energy densities of hessence and background matter are given by 

n h = x 2 +y 2 ~v 2 , n m = z 2 , (29) 
respectively. The EoS of hessence and the effective EoS of the whole system are 



Ph X 2 -V 2 -y 2 Ph+Pm 2 2 2 2 

w h = — = — — — w eff = ■ =x -v -y + 7 - 1 )z , 30 

p h x z - v z + y l p h + pn 



respectively. One can see that t«ft>-lasi 2 >w 2 , while wt < — 1 as x 2 < v 2 . Finally, it is worth noting 
that y > and z > by definition, and in what follows, we only consider the case of expanding universe 
with H > 0. 

It is easy to see that Eqs. (|21ll - H25|l become an autonomous system when the potential V(4>) is chosen 
to be an exponential or (inverse) power law potential and the interaction term C is chosen to be a suitable 
form. Indeed, we will consider the model with an exponential or (inverse) power law potential in Sec. IIVI 
and Sec. respectively. In each model with different potential, we consider four cases with different 
interaction forms between hessence and background perfect fluid. The first case is the one without 
interaction, i.e. C — 0. The other three cases are taken as the most familiar interaction terms extensively 
considered in the literature: 

Case (I) C = 0, 

Case (II) C — aKp m cl), 

Case (III) C = WHptot = 3(3H [p h + p m ) , 

Case (IV) C = 3r/Hp m , 

where a, (3 and r\ are dimensionless constants. The interaction form Case (II) arises from, for instance, 
string theory or scalar-tensor theory (including Brans-Dicke theory) [33l Isfij . The interaction forms 
Case (III) [2a| and Case (IV) [22| are phenomenally proposed to alleviate the coincidence problem. 

In the next two sections, we first obtain the critical points (x,y,z,u,v) of the autonomous system by 
imposing the conditions x' = y' = z' = vl = v' = 0. Of course, they are subject to the Friedmann 
constraint, i.e. x 2 + y 2 + z 2 — v 2 = 1 . We then discuss the existence and stability of these critical points. 
An attractor is one of the stable critical points of the autonomous system. 



G 



IV. 



MODEL WITH EXPONENTIAL POTENTIAL 



In this section, we consider the hessence model with an exponential potential 



(31) 



where A is a dimensionless constant. Without loss of generality, we choose A to be positive, since we can 
make it positive through field redefinition <f> — ► —(f) if A is negative. In this case, Eqs. (|21Jl - l|25|l become 

3 



3x (^x 2 


-v 2 - 






-v 2 - 


-lz 2 
2 


3z [x 2 




-lz 2 
2 


—xu 2 , 






3v (x 2 


-« 2 H 


-'¥ 



1 



+ C 2l 



(32) 

(33) 

(34) 
(35) 
(36) 

where C\ and C 2 are defined by Eq. (|26l) and depend on the interaction form C. In the following 
subsections, we will consider different interaction forms C mentioned in the end of Sec. IIIII 

To study the stability of the critical points of Eqs. (|32|I - H36|I . we substitute linear perturbations x — ► 
x, + 5x, y — > y + Sy, z — > z + 5z, u — > u + 5u, and v — > v + Sv about the critical point (x, y, z, u, v) into 
Eqs. (|32() - (|36|l and linearize them. Note that these critical points must satisfy the Friedmann constraint, 
y > 0, z > and requirement of x, y, z, u, v all being real. Because of the Friedmann constraint l)28[l. 
there are only four independent evolution equations: 



Sx' 
Sy' 

5Z' : 

Su' 



3y 2 



3 1 



D z 2 + 2ux 



Sx 



+ ( V6A - 6x 



2uj y5y 



+ [3 (7 - 2) x - 2u] zSz - (x 2 + y 2 + z 2 -l)5u- Sd, 



2 

-QyzSy 



-u 2 Sx 



Xy5x 



1-W 

_ 7 
2 



6y + 3{~f-2)yz5z, 



6z + 6C 2 , 



2xuSu, 



(37) 
(38) 

(39) 
(40) 



where 5C\ and bCi are the linear perturbations coming from G\ and C 2 , respectively. The four eigenvalues 
of the coefficient matrix of the above equations determine the stability of the critical point. 



A. Case (I) C = 



C = means no interaction between hessence and background matter. In this case, one has C\ = C% = 
0. It is easy to find out all critical points (x, y, z, u, v) of the autonomous system l|32ll -H36 |l . They are 
required to be real and satisfy the Friedmann constraint and y > 0, z > 0. We present them in Table 
Next we consider the stabilities of these critical points. In this case, 8C\ = 6C 2 = 0. Substituting SCi, 
SC 2 and the critical point (x,y,z,u,v) into Eqs. i|3T |l -l|l0" |l . we find that Points (E.I.I), (E.I.2), (E.I.3) 
and (E.I. 4) are always unstable, while Point (E.I. 5) and Point (E.I. 6) exist and are stable under condition 
A > \fWy and A < a/37, respectively. 

The late time attractor (E.I. 5) has 

37 37 

Vl h = —, fl m = l-— , 10/, =-1+7, W e ff = -1+7, (41) 

which is a scaling solution. The late time attractor (E.I. 6) has 

A 2 A 2 
Sl h = l, £l m = 0, w h = -l + —, w eff = -! + —, (42) 



which is a hessence-dominated solution. Note that their EoS are all larger than — 1. 
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Label 


Critical Point (a;, y, z, u, v) 


E.I.I 


> 1, 0, 0, 0, ±Vx 2 - 1 


E.I. 2 


0, 0, 1, any, 


E.I.3 


±1, 0, 0, 0, 


E.I.4 


^,0, 0, 0,±y/&-l 


E.I.5 




E.I.6 


^, ^/l-f ,0, 0,0 



TABLE I: Critical points for Case (I) C = in the model with exponential potential. 



Case (II) C = anpr, 



In this case, C\ 



az 2 and Co 



-axz. The physically reasonable critical points of the au- 



tonomous system H32|+()36|) are summarized in Table ITT1 Next let us consider the stability of these critical 

ax5z. Substituting 8C±, 8C2 and the 



^azSx 



points. In this case, 8C\ — \/6azSz and 8C1 
critical point (x,y,z,u,v) into Eqs. H37|l - (|40[l . we find that Point (E.II.l) exists and is stable under 
condition a < and x > max|l, ^/f-Tpj; Point (E.II.2p) exists and is stable under condition 

ot < \l 55(7 — 2) and A > ^/6; Point (E.II.2m) is always unstable; Point (E.II.3) exists and is stable under 



2 

condition A > 1/67(2 - 7) and \/§(7-2) < a < </§(2-7) and 



-\- y /\ 2 -6~ t (2-~ t ) -A+^/A 2 -6 7 (2- 7 ) . 

2 ^ 2 ' 



Point (E.II.4) exists and is stable under condition A < -s/6 and a < — l) A; Point (E.II.5) exists and 

is stable under condition > a > max |^/|(7 — 2), — l) a|; Point (E.II.6) exists and is stable in 

a suitable parameter-space 38]; and Point (E.II.7) exists and is stable under condition A < \/6 and 
a < -?■ — A. The corresponding fi/j, ^m, Wh and w e ff of the attractors are presented in Table ITT1 as well . 
Again, we see that their EoS of hessence and effective EoS are always larger than —1 [note that in Point 
(E.II.6), a(X + a) + 37 > and a + A > arc required by the existence of its corresponding y and z, 
respectively]. This implies that the big rip singularity will not appear in this case. 



Label 



Critical Point (a;, y, z, u, v) 



U'h 



"eff 



E.II.l 
E.II.2p 
E.II.2m 
E.II.3 
E.II.4 
E.II.5 
E.II.6 
E.II.7 



3 -1 

2 \ + a ' 



x 2 > 1, 0, 0, 0, ±Vx 2 - 1 
+1, 0, 0, 0, 
-1, 0, 0, 0, 



2 a n /1 2a 2 r. r, 

3 T -2' U > Y 3(7 _2)2 ' U > U 

V6 



0, 0, 0, ± 



--^ o,o,o, ±v /^ 



2a 2 -3(7-2)7+2aA 



2(<* + A)^ 



A 

V6 



y/l-%,0, 0,0 



1 
1 

2c? 



3(7-2) 2 



2c 2 
3( 7 -2) a 



a{\ + a) + 3"l A(A + q)-37 
(A + ap (\ + a) 2 



-1 + 



37" 



ck(A + ck) + 37 



-1 



A^ 

3 



-1 + 7 + 



-1 + 
-17 



2a 2 



3(2-7) 



TABLE II: Critical points for Case (II) C = a.Kp m 4> in the model with exponential potential. 
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C. Case (III) C = 3(3Hp tot = 3(3H {p h + p m ) 

In this case, C\ — and C2 = The physically reasonable critical points of the autonomous 

system (|32|) - (j36[) are presented in Table ITTT1 where 



4(3 
2-7' 



(43) 



Points (E.III.4), (E.III.5) and (E.III.6) are the three real solutions of 




x + 1, 



/?, (44) 

(45) 
(46) 



with u — v — and requiring y > 0, z > by definition. Although the real roots of the Eqs. (|44H 1461) 
are easy to obtain, we do not present them here, since they are long in length. Note that Point (E.III.3) 
exists only when > and 7 > 2, which is beyond the range < 7 < 2. 

Next, we consider the stabilities of these critical points. In this case, 8C\ = —^(3x~ 2 5x and 6C2 = 
— \(5z^ 2 5z. Substituting SCi, 8C2 and the critical point (x,y,z,u,v) into Eqs. l|3*7|l - (|4T1| . we see that 
Points (E.III.lp), (E.III.lm) and (E. III. 2m) are always unstable, while Point (E.III.2p) exists and is stable 
under condition > (3 > (7 - 2)/4 and A > [12 + 3(1 + r 3 )(-2 + 7)] /J2y/5y/l - r 3 ). Points (E.III.4), 
(E.III.5) and (E.III.6) exist and are stable in proper parameter-space [3H| - 

The late time attractor (E.III.2p) has 

ft h = -{l-r 3 ), O m = i(l + r 3 ), w h = l, w eff = 1 + i (7 - 2) (1 + r 3 ) , (47) 
which is a scaling solution. The late time attractors (E.III.4), (E.III.5) and (E.III.6) have 



y 



w h 



1 



2x 2 



W eff 



1 



Ax. 



(48) 



They are all scaling solutions. It is easy to see that their EoS of hessence and effective EoS are all larger 
than —1 [note that < 7 < 2 and < r 3 < 1 for Point (E.III.2p), while x > is required by Eq. for 
Points (E.III.4), (E.III.5) and (E.III.6)]. 



Label 



E.III.lp 
E.III.lm 
E.III.2p 
E.III.2m 
E.III.3 



Critical Point (x, y, z, u, v) 



[i(l+r 3 )] 1/2 ,0, [i(l-r 3 )] 1/2 ,0,0 
-[i(l + r 3 )] 1/2 ,0, [i(l-r 3 )] 1/2 ,0,0 

[i(l-r 3 )] 1/2 ,0, [i(l + r 3 )] 1/2 ,0,0 
-[i(l-r 3 )] 1/2 ,0, [I(l + r 3 )] 1/2 ,0,0 
0, ± 



A ' 



-2+7' 



-! + ■ 



M 1 6 
4+27 A 7 



Note: For Points (E.III.4), (E.III.5) and (E.III.6), see text 



TABLE III: Critical points for Case (III) C = 3/3Hp to t = 3/3H (ph + pm) in the model with exponential potential. 
r 3 is given in Eq. I|48^ . 
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D. Case (IV) C = Z-qHp m 

In this case, C\ — ^rjx~ 1 z 2 and C2 = \ i)Z. We show the physically reasonable critical points of the 
autonomous system l|32[) (1361) in Tabic II VI where 



-3 7 3 + 6 7 2 (1 + t?) - 3 7? ?(4 + 77) + 2^(37? + A 2 ) 

2 7 A 2 ' l4yj 



_ / ( 7 - 7?)A 2-3( 7 -7 ? ) 2 



(50) 



To study the stability of these critical points, we obtain 8C\ 



\rjx 2 z 2 Sx + 3r]x 1 z5z and 8C2 = %r\bz 



by linearizing Ci and C2. Substituting <5Ci, <5C2 and the critical point (x,y,z,u,v) into Eqs. 137|l (|40|1 . 
wc find that Point (E.IV.l) exists and is stable under condition 77 < —2 + 7 and x > max |l, j; Point 



(E.IV.2p) exists and is stable under condition > r\ > 7 - 2 and A > y 5(7 — ry) y 2 ^- 2 -; Point (E.IV.2m) 

is always unstable; Point (E.IV.3p) exists and is stable under condition 7/ < 7 — 2 and A > v6; Point 
(E.IV.3m) is always unstable; Point (E.IV.4) exists and is stable under condition 7/ < 7 — 2 and A < \/6; 

Point (E.IV.5) exists and is stable under condition 77 < 7 and A < min|\/6, \/3( 7 — ?/) j; and Point 

(E.IV.6) exists and is stable in proper parameter-space |3S| . 

The corresponding f2^, f2 TO , tw^ and w e // of the attractors are presented in Tabic llVl as well. Once 
again, we see that their EoS of hessence and effective EoS are all larger than —1 [note that 7/ < is 
required by the existence of its corresponding x for Point (E.IV.2p), while 7 — 77 > is required by 
Eq. (jHOJ for Point (E.IV.6)]. 



Label 



Critical Point (x, y, z, u, v) 



77 



Weff 



E.IV.l 

E.IV.2p 
E.IV.2m 
E.IV.3p 
E.IV.3m 
E.IV.4 
E.IV.5 
E.IV.6 



> 1, 0, 0, 0, ±Vx 2 - 1 



7-2 ' 



0, 



7-2 ' u ' 



1-^5,0,0 



+1, 0, 0, 0, 
-1, 0, 0, 0, 



^,0, 0, o,±j£ 



3 i-v r r n n 



n 1 _ v 

7-2 7-2 



-1 + — 

1 ~ 3 



-1 + 



3(7~'7) 2 
A 2 (l-rj) 



-1 + 7 — 77 



-1 + V 



-1+7 — 77 



TABLE IV: Critical points for Case (IV) C = 3r]Hp m in the model with exponential potential. r y and r z are 
given in Eqs. I|49|l and 1500 , respectively. 



V. MODEL WITH (INVERSE) POWER LAW POTENTIAL 

In this section, we consider the hessence model with a (inverse) power law potential 

V{<j>) = V W , (51) 
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where n is a dimensionless constant. V(<f>) is a power law potential when n > while it is an inverse 
power law potential when n < 0. In this case, Eqs. H21|) - (|25[) become 



' o ( 2 2 

x = 3x \x — v 
V = 3y (x 2 



7 2 



z — 1 ) — uv — —uy — Ci, 



I o I 2 2 

Z = 6Z \x — V 

.2 



V 



-xu 



iv (x 2 



V 



2 2 



2 



Co 



SOT, 



(52) 
(53) 
(54) 
(55) 
(56) 



where C\ and C2 are dehned in Eq. (|26|) and depend on interaction form C . In the following subsections, 
we will consider four interaction forms C given in the end of Sec. II I II 

To study the stability of the critical points of Eqs. (|52|I - H56|) . we substitute linear perturbations x — > 
x + Sx, y —>■ y + Sy, z — * z + 5z, u — + u + 5u, and v — > v + <5w about the critical point (x, y, z, u, S) into 
Eqs. (|52|) - (|56|l and linearize them. Note that these critical points must satisfy the Friedmann constraint, 
y > 0, z > and requirement of x, y, z, u, v all being real. Because of the Friedmann constraint (|28|l . 
there are only four independent evolution equations: 



5x' 



-1) : 
1 



2mx| Sx — [6xy + (2 + n) uy] Sy + [3 (7 — 2) x — 2u] zSz 



Su — 5C\, 



—yuSx + 3 



Sy' 
Sz' 

Su 1 = —u 2 Sx — 2xu5 



6yz5y + | (7 - 2) . 



1-1 
2 



n __ 
—ux 
6 



$y + 3 (7 — 2) yzbz + —xySu, 



i-D-f+G-i)*]}**'*, 



(57) 

(58) 

(59) 
(60) 



where 5C\ and 8C2 are the linear perturbations coming from C\ and C2, respectively. Again, the four 
eigenvalues of the coefficient matrix of the above equations determine the stability of the critical points. 



A. Case (I) C = 

In this case, C\ — C2 — 0, namely there is no interaction between hessence and background matter. It 
is easy to find out all critical points (x, y, z, u, v) of the autonomous system l(^ )l -l(5^ )l . We present them 
in Table|V| Checking the eigenvalues of Eas. ^ -fffl ty . we find that Points (P.I.I), (P.I.2), and (P.I.4) are 
always unstable, while Point (P.I. 3) is always stable. The unique late time attractor (P.I. 3) has 

Q/, = 1, O m = 0, w h = -1, w eff = -1, (61) 

which is a hessence-dominated solution and it is an asymptotical de Sitter attractor. 



Label 


Critical Point (x, y, z, u, v) 


P.I.I 


x z > 1, 0, 0, 0, ±yjx 2 - 1 


P.I.2 


0, 0, 1, any, 


P.I.3 


0, 1, 0, 0, 


P.I.4 


±1, 0, 0, 0, 



TABLE V: Critical points for Case (I) C — in the model with (inverse) power law potential. 
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Case (II) C = anp-n 



In this case, C\ — y i|<X2 and C2 — y ^axz. We summarize the physically reasonable critical points 

of the autonomous system (I52(l - (|56[1 in Table IvTl Next, we consider the stabilities of these critical points. 

In this case, 8C\ = \/6azSz and SC2 = y^azSx + y^axSz. The four eigenvalues of the coefficient 

matrix of Eqs. dHZI dHOl) tell us that Points (P.II.l), (P.II.2), (P.II.3) and (P.II.4) are always unstable, 
while Point (P.II.5) is always stable. The unique late time attractor (P.II.5) has 



fl h = 1, Q m = 0, w h = -1, W ef f = -1, 
which is a hessence-dominated solution and it is an asymptotical de Sitter attractor. 



(62) 



Label 



Critical Point (x, y, z, u, v) 



P.II.l 
P.II.2 
P.II.3 
P.II.4 
P.II.5 



> 1, 0, 0, 0, ±Vi 



±1, 0, 0, 0, 




2 _a 

3 7-2 ' v '' 

± 

2 a ' ' ' ' 



2q 2 



3(7-2) 2 ' 



0, 



1 _l_ 3(7~2) 2 

_I H 9^— 



TABLE VI: Critical points for Case (II) C = aKp m <f> in the model with (inverse) power law potential. 



C. Case (III) C = 3f3Hp tot = 3/3 H (p h + p m ) 

In this case, C\ = and C2 = The physically reasonable critical points of the autonomous 

system (IK2"l) - (f5ri)) are summarized in Table IVlTl There, r 3 is defined by Eq. (|4*3|l . Substituting 8C\ — 
— ^(3x~ 2 5x, SC2 = — ^{3z~ 2 5z, and the critical point (x, y, z, u, v) into Eqs. (|57|l l|6()(l . we find that no 
stable attractor exists in this case. 



Label 



Critical Point (x, y, z, u, v) 



PIII.lp 
P.III.lm 
P.III.2p 
P.III.2m 



[i(l + r 3 )] 1/2 ,0, [i(l-r 3 )] 1/2 ,0,0 
-[I(l + r 3 )] 1/2 ,0, [I(l-r 3 )] 1/2 ,0,0 
[I(l^r 3 )] 1/2 ,0, [i(l+r 3 )] 1/2 ,0,0 



-[i(l-r 3 )] 1/2 ,0, [i(l + r 3 )] i/2 ,0,0 



1/2 



TABLE VII: Critical points for Case (III) C = Z(3Hptot = 3(3H (p h 
potential. r 3 is given in Eq. 143H . 



pro) in the model with (inverse) power law 



D. Case (IV) C = 3nHp„ 



In this case, C\ = \i}x 1 z 2 and C2 = §?7^, we have the physically reasonable critical points shown in 
Table IVlHl Substituting 8C\ — —^rjx~ 2 z 2 Sx + 3i]x~ 1 z8z, 8C2 — ^Sz, and the critical point (x, y, z, u, v) 
into Eqs. H57|) - H60[) . the four eigenvalues of the coefficient matrix of the resulting equations tell us that 
no stable attractor exists as well in this case. 
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Label 


Critical Point (a;, y, z, u, v) 


P.IV.l 


> 1, 0, 0, 0, ±Vx 2 - 1 


P.IV.2 


±V^. o, ^1-^,0,0 


P.IV.3 


±1, 0, 0, 0, 



TABLE VIII: Critical points for Case (IV) C = SrjHp m in the model with (inverse) power law potential. 

VI. NO BIG RIP IN HESSENCE MODEL 

Obviously, we can see from Sec. IIVI and Sec. that Wh and w e ff of all stable late time attractors 
are larger than or equal to —1 for these two particular models with exponential or (inverse) power law 
potential and the interaction term C is chosen to be four different forms. No phantom-like late time 
attractor with Wh or w e // less than —1 can exist in the hessence model of dark energy. However, one 
may wonder whether this observation depends on the forms of potential V((j>) and the interaction term 
C. Therefore, it is interesting to investigate this feature in a general case. 

Let us come back to the most general equations system, i.e. Eqs. H21|l - (|25|) . Now, we leave V(4>) and 
C undetermined, except for assuming that they can make the equations system closed, in other words, 
Eqs. I|21|) - (|25|l is an autonomous system. The critical point (x, y, z, u, v) of the autonomous system 
satisfies x' — y' — z' — u' — v' — 0. Of course, they are also subject to the Friedmann constraint, 
i.e. x 2 + y 2 + z 2 — v 2 = 1. From Eq. J21J, u' — 0, we have xu — 0. Substituting this into Eq. (J2EJl, 
v' = 0, we have either v — or x 2 — v 2 + 7;z 2 — 1 = 0. If the latter holds, from Eq. (|29|l . we have 
Q m = z 2 — — [l — (i 2 — v 2 )] . Physics requires £l m < 1. Therefore x 2 — v 2 > is required for the case 
< 7 < 2. We see from Eq. l|!?0)l that Wh > —1 in this case. If x 2 — v 2 + ^z 2 — 1 ^ 0, we then have v — 
as mentioned above. Obviously, from Eq. (|3()|) . u>h > — 1 is inevitable. On the other hand, we can see no 
big rip in the hessence model as follows. Since x 2 — v 2 > for both cases mentioned above, we can see 
from Eq. (|27|l that H < 0, which implies w e ff > —1, cf. Eq. I]16p. Therefore, we conclude that Wh > — 1 
and w e ff > — 1 always hold for all critical points of the autonomous system l|21|) - (|25|l . No phantom- like 
late time attractor with EoS less than —1 can exist. This result is independent of the form of potential 
V((j)) and interaction form C. 

Therefore, EoS less than —1 is transient in the hessence model. Eventually, it will go to quintessence- 
like attractors whose EoS is larger than —1 or asymptotically to de Sitter attractor whose EoS is a 
constant —1. Thus, the big rip will not appear in the hessence model. 



VII. CONCLUSION AND DISCUSSIONS 



In this work, the cosmological evolution of hessence dark energy is investigated. We considered two 
models with exponential and (inverse) power law potentials of hessence respectively, and investigated the 
dynamical system for the four cases with different interactions between hessence and background perfect 
fluid. By the phase-space analysis, we find that some stable attractors can exist, which are either scaling 
solutions or hessence-dominated solutions with EoS larger than or equal to —1. No phantom- like late 
time attractors with EoS less than —1 can exist. We have shown that this essential result still holds in a 
general case beyond particular potentials and interaction forms. Thus, the big rip will not appear in the 
hessence model. 

Our result of the hessence model is very different from that of the quintom model studied in [22, H3| , 
where the phantom-dominated solution is the unique late time attractor and the big rip is inevitable. 
The difference between our result and that of [2^, |2^| is not due to the interaction between the hessence 
and background matter, since our conclusion still holds for the case C = 0, that is, the case without 
the interaction. The difference should be due to the form of the interaction between <f>i and 02 [cf. 
Eqs. 0-0]. Guo et al and Zhang et al [2j| only studied the cosmological evolution of the quintom 
model without direct coupling between 4>\ and 4>2, i-e. V(4>i,(j)2) = V ( f >1 + V<f, 2 = V < f >1 oenp(—X ( f >1 K,4 > i) + 
V0 2 o exp(— \<p 2 k4>2), and with a special interaction between 0i and 02, i-e. V(tf>i, 4>2) — V^ t + V^ 2 + Vint 
and Vi n t ~ (V ( j >1 Vfa ) 1 ^ 2 , respectively. As pointed out in 28], in the hessence model, the potential is 
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imposed to be the form of V(4>), or equivalently, V{4>\ — 4>\) in terms of <f>\ and </>2 [cf. Eqs. Q-Q]. 
Except for the very special case with V{4>) ~ <fr 2 , the two fields 4>i and 4>2 are coupled in the general 
case, and the interaction between (f>\ and <p2 is quite complicated, rather than the very special interaction 
considered in [^J. 

Another issue is about the fate of our universe. From our result, the present super-acceleration of 
universe, i.e. H > 0, is transient and our universe can avoid the fate of big rip. It can either accelerate 
forever (iw e // < —1/3) or come back to deceleration (w e ff > —1/3), which is determined by the model 
parameters. 
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